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We examine simulations of the formation of domain walls, cosmic strings, and monopoles on a 
bJQ. \ 

cubic lattice, in which the topological defects are assumed to lie at the zeros of a piecewise constant 

1, 2, or 3 component Gaussian random field, respectively. We derive analytic expressions for the 

corresponding topological defect densities in the continuum limit and show that they fail to agree 

with simulation results, even when the fields are smoothed on small scales to eliminate lattice 

effects. We demonstrate that this discrepancy, which is related to a classic geometric fallacy, is 

QQ . due to the anisotropy of the cubic lattice, which cannot be eliminated by smoothing. This problem 

can be resolved by linearly interpolating the field values on the lattice, which gives results in good 

■ agreement with the continuum predictions. We use this procedure to obtain a lattice-free estimate 

(for Gaussian smoothing) of the fraction of the total length of string in the form of infinite strings: 



faa = 0.716 ± 0.015. 



On 

pu,: 

43 I. INTRODUCTION 



Topological defects arise when the manifold of the degenerate vacuum states of a field is not simply connected. 
(See Ref. Q] for a review). In the simplest set of models, domain walls arise from a one-component field with two 
■ degenerate vacuum states, cosmic strings correspond to a two-component field in which the degerate vacua form a 
circle in the field space, and monopoles arise from a three-component field for which the degenerate vacua lie on the 
surface of a sphere. 

The formation of topological defects can be simulated by laying down the appropriate field on a cubic lattice, and 
identifying the corresponding topological defects with zeros of the field. Although the most effort has gone toward 
simulations of cosmic string formation 01 - |3L similar investigations have also been undertaken for the study of 
domain walls |[(| - JlJ] and monopoles jl3| - [14(|. 

In this paper, we point out a potential problem for some such lattice-based simulations. Even when the fields which 
give rise to the defects are smoothed on a scale larger than the lattice spacing (as in Ref. ||), residual lattice effects 
remain. These are due to the fact that a continuous surface cannot be distorted to lie on the edges of a lattice without 
also distorting the area of the surface. 

In the next section, we discuss our topological defect simulations, and we provide analytic expressions for the 
density of walls, strings, and monopoles. These analytic predictions are compared with the results of our numerical 
simulations, and they are found to disagree. In Sec. 3, a second set of analytic predictions, which does agree with the 
simulations, is presented, and we explain the origin of the discrepancy, which is closely related to a classic geometric 
fallacy. A second set of simulations using linear interpolation of the fields is presented, and these are shown to agree 
with analytic continuum predictions. The implications are discussed in Sec. 4. Our interpolation scheme produces a 
lattice-free estimate for the fraction of string length in infinite strings (with Gaussian smoothing): foo — 0.716±0.015. 
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II. THE PROBLEM 



In an earlier paper [gj we simulated the formation of cosmic strings from correlated fields on a cubic lattice, 
smoothing the field on small scales to eliminate lattice effects. Here we extend these simulations to domain walls and 
monopoles. 

Consider first the case of domain walls. We can assign a real scalar field <p to each of the cells of a cubic lattice, 
and the location of the domain wall is then identified with faces of the lattice across which the value of the field 
crosses zero. To simulate cosmic strings rather than domain walls, we use a complex scalar field, or equivalently, a 
two-component real field <f>. Then the zeros of <fi will lie on the edges of the lattice, which we then identify as the 
location of the strings. Finally, to simulate monopoles, we take <fi to be a three-component real field; the zeros of <p 
then form a set of disconnected points lying on the vertices of the lattice; this gives the location of the corresponding 
monopoles. 

Following the procedure in Ref. we take to be a Gaussian random field. Then we can exploit the fact that 
the components of an N— component Gaussian field, <p\, 4>2, ■■■4>n, are themselves independent real Gaussian fields. 
Hence, our model for cosmic strings is equivalent to using two independent real Gaussian fields (which can be taken 
to be </>i and </>2). The zeros of these two fields form two independent sets of surfaces, and the strings lie at the 
intersection of these surfaces. Similarly, our field configuration for monopoles can be treated as three independent 
real Gaussian fields, <f>i, 02, and </>3. The zeros of these three fields form three sets of surfaces, and their intersection 
is a set of disconnected points, giving the location of the monopoles. 

A Gaussian field is completely determined by its power spectrum P(k), defined by 

P(k) = J d 3 r e 4k - r (0(x) • 0(x + r)). (1) 

Following the procedure in Ref. (^] , we assume a power-law power spectrum of the form 

P{k) oc k n . (2) 

Note that this is equivalent to taking each of the components (j>i to have a power spectrum given by equation (||) . We 
add a backround non-zero mean field to the box to simulate the effect of long-wavelength modes (see Ref. || for the 
details) . We also smooth the field 4> on small scales with a Gaussian window function 

W(r) =exp(-r 2 /2r2), (3) 

where the smoothed field s (x) is the convolution of 0(x) with VF(r): 

s (x) = J d 3 r 0(x + r)VF(r). (4) 

This Gaussian smoothing changes the power spectrum from P{k) to P s (k), given by 

P.(k) = P(k)e~ k2r Z (5) 

(Again, equation (|]) is equivalent to smoothing each of the components (f>i with W(r), and equation (||) can similarly 
be applied separately to the power spectrum for each component.) The effect of smoothing is to reduce the magnitude 
of the small-scale fluctuations by averaging them out over the window function. One might hope that if ro is taken 
to be larger than the lattice spacing, lattice effects can be reduced or eliminated entirely We will see that this is not 
the case. 

Given this model, it is possible to calculate the topological defect density analytically. Our starting point is a set 
of results due to Ryden jl5| concerning the properties of level-crossing surfaces in Gaussian random fields. Consider 
a real Gaussian field with arbitrary power spectrum P(k) and rms fluctuation a. The regions of space for which the 
field has the value va form a set of two-dimensional surfaces. Ryden jfl| defined the quantity Ns(v) to be the mean 
area per unit volume of these surfaces. The intersection of this set of surfaces with an arbitrary plane produces a set 
of curves, with mean length Ni(y) per unit area of the plane. Finally, we can consider the intersection of this set of 
surfaces with an arbitrary straight line; the quantity N\{v) is the number of intersections of this line (per unit length 
of the line) with the set of surfaces. Ryden showed that Il5[ 

N 3 (y) = -N 2 (u) = 22V x (i/), (6) 

7T 
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where 
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Nl (,) = -^=(kr /2 e-^/ 2 , (7) 



and the value of (fc 2 ) depends on the power spectrum as 

_ JP(k)k*dk 

{ } ~ jP(k)k 2 dk W 

These results can be used to derive analytic expressions for the topological defect density in our model. Consider 
first the area per unit volume, (A/V), for domain walls. In terms of the Ryden notation, it is obvious that 

(A/V)=N 3 (Q) = ^ 7 ={k 2 ) 1/2 (9) 

7TVo 

(Although we confine our attention to unbiased field configurations, i.e., configurations for which positive and negative 
values of 4>i are equally likely, so that v = 0, our results can easily be generalized to biased configurations by taking a 
nonzero value for v). Now consider the case of cosmic strings. The strings can be considered to lie at the intersection 
of two sets of v — surfaces, corresponding to <f>\ = and 4>2 = 0, respectively. Suppose that we are sitting on a 
surface corresponding to <\>\ = 0, with area per unit volume A^O), and we wish to calculate the length per unit area 
of the intersection of this surface with the (j>2 = surface. Although the </>\ — surface is highly irregular, if we take 
a small enough patch, it will be locally flat (assuming the field is smoothed on small scales), so that the length per 
area of this intersection is just N^ifS). Hence, the total length per unit area, (L/V), of cosmic string is just 

(L/V) = N 2 (0)N 3 (0) = ^(fc 2 ) (10) 

Note that this result agrees exactly with the analytic expression for (L/V) derived by Vishniac et al. using 
different methods. Finally, we consider the case of monopoles, which are taken to lie at the intersection of the surfaces 
corresponding to the zeros of three independent Gaussian fields. The intersection of two of these sets of surfaces 
(4>\ = and 02 = 0) defines a set of curves with length per unit volume given by equation (|l^). While these curves 
are quite irregular, they are locally straight on small enough scales (assuming the fields are smoothed on small scales), 
so that their intersection with the third set of surfaces (4> 3 = 0) gives Ni (0) as the number of intersections per unit 
length of the curves. Then the monopole density per unit volume, (N/V), is 

(N/V) = A 1 (0)A 2 (0)A 3 (0) = ^^(fc 2 ) 372 (11) 

Although our results are applicable to arbitrary power spectra, for definiteness we now restrict our attention to the 
case n = 0, which corresponds to uncorrelated fields. This is the most physically relevant case, because a causal phase 
transition will lead to fields which are uncorrelated on scales larger than the horizon, while the smoothing length can 
be identified with the scale over which the field tends to be correlated. For n = with Gaussian smoothing, equations 
(§),(§, and © yield 

< fc2 > = ^ (12) 



Substituting this result into equations (^|), ( |l0| ) and (|ll|), we obtain 

(A/V) = — (walls), 
7rr 

{L/V) = 2^1 (brings), 

WV) = ^7^3 ( mon °P° les )- ( 13 ) 

We now compare these results to numerical simulations. Using the procedure outlined above, we have simulated the 
formation of domain walls, cosmic strings, and monopoles on a 128 3 lattice, using 1, 2, and 3 independent Gaussian 
random fields, respectively. For the case of cosmic strings, we use two different cubic lattices, staggered with respect 
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to each other by 1/2 of a lattice spacing in the x, y, and z directions (i.e., the vertices of one lattice lie at the centers 
of the cells of the other lattice). We then set down independent Gaussian fields on each lattice, and the strings are 
taken to lie at the intersection of the zeros of the two fields. For monopoles, we take three different cubic lattices, 
staggered with respect to each other by 1/3 of a lattice spacing in the x, y, and z directions, and the monopoles are 
taken to lie at the intersection of the three surfaces defined by </>i = 0, (f>2 — 0, and $3 = 0. 

In each case, we produce 4 different realizations to obtain a mean defect number density and standard deviation. 
In Fig. 1, we show the domain wall area per unit volume, (A/V), in our simulations as a function of smoothing length 
ro. In Fig. 2 we give the string length per unit volume, (L/V), as a function of ro, and in Fig. 3 we give the number 
of monopoles per unit volume, (N/V), as a function of ro. (Here we make no distinction between monopoles and 
antimonopoles, although we will do so later). Also in these figures, we show as solid lines the theoretical values for 
(A/V), (L/V), and (N/V) from equations (|l3|). The disagreement is obvious. Clearly the lattice-based simulations 
fail to produce the true continuum number densities for topological defects. 
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FIG. 1. The total domain wall area per unit volume, (A/V), multiplied by the smoothing length ro, as a function of ro, 
for Gaussian smoothing, where ro is measured in units of the lattice spacing. The solid line gives the analytic prediction for 
(A/V)ro for this model in the continuum limit. The dotted line is the lattice-based analytic prediction. 
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FIG. 2. The total string length per unit volume, (L/V), multiplied by the square of the smoothing length ro, as a function 
of ro, for Gaussian smoothing, where ro is measured in units of the lattice spacing. The solid line gives the analytic prediction 
for (L/V)ro for this model in the continuum limit. The dotted line is the lattice-based analytic prediction. 
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FIG. 3. The total number of monopoles per unit volume, (N/V), multiplied by the cube of the smoothing length ro, as a 
function of ro, for Gaussian smoothing, where ro is measured in units of the lattice spacing. The solid line gives the analytic 
prediction for (N/V)rQ for this model in the continuum limit. The dotted line is the lattice-based analytic prediction. 
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III. THE SOLUTION 



Although our analytic expressions for the topological defect densities given in the previous section fail to agree with 
the results of our numerical simulations, we can use a different set of arguments to derive "lattice-based" analytic 
expressions which do agree with the lattice simulations. Consider first the case of domain walls, and consider a single 
face between two cells on the lattice. This face will be occupied by a domain wall if the cell on one side of it has 
<t> > and the other cell has < 0, or vice versa. This occupation probability, p cc, can be derived from Sheppard's 
theorem [fL7| when is a Gaussian field. We obtain 



1 -i 

Voce = - COS 



(14) 



where r is the spacing between the lattice cells, and £ is the two-point correlation function, which is the Fourier 
transform of the power spectrum. For an n — power spectrum smoothed with our Gaussian window function, we 
have £(r)/£(0) = cxp(-r 2 /4rg), and 



1 -i 

Voce = — COS 



-r 2 /4r= 



(15) 



For rp > r (i.e., when the smoothing length is large compared to the lattice spacing) this reduces to 



Voce = — t=— (16) 

This result can be used to derive not only the domain wall density, but the string and monopole densities as well. 
For the case of domain walls, each cell is bounded by six faces, each of which is shared between two cells, so that 
the total area per unit volume is (A/V) — 3p occ . In our string simulation, a string forms at the intersection of two 
different lattice faces, where the two lattices are staggered by half a lattice spacing. For a single cell on one lattice, 
there are 24 such possible intersections, each with length r/2 (where r is the lattice spacing) and each shared by two 
cells. Hence, the total length per unit volume is (L/V) = 6p„ cc . Finally, for the monopoles, a monopole forms at 
the intersection of three lattice faces on three different staggered lattices. On a single cell of one lattice, there are 12 
possible intersection sites, each of which is shared between two cells. Thus, (N/V) = Qvlcc- 

Combining these results with the value for p occ given in equation (|l6|), we obtain the following values for the 
topological defect densities: 

(A/V) lat = — (walls), 
7r\/2ro 
3 

(L/V)i at = (strings), 



(N/V)i at = 3 3 (monopoles), (17) 

7T \J LVr. 



where the lot subscript denotes the fact that this analytic calculation is lattice-based, rather than derived from a 
continuum calculation. 

These values for the defect densities are shown in Figs. 1-3 as dotted lines. As r becomes sufficiently large 
compared to the lattice spacing, the numerical results converge to a defect density in good agreement with equations 
(|i"7|), although the fluctuations between realizations also increase with tq. Clearly, our results in equations ([17]) give 
the correct defect densities to compare with our lattice-based simulations. 

However, this leaves an embarrassing question: why do our simulations fail to agree with the continuum calculations 
given in equations ( |l3| ) and in reference fill ? Even without reference to the numerical simulations, we have two 
predictions for the defect densities, i.e. equations (|T^ ) and (|l7|), which disagree with each other. What is the origin 
and significance of this discrepancy? 

In fact, the problem here is related to a classic geometric fallacy flq| . Consider an isosceles right triangle with 
sides of unit length. The hypotenuse can be represented as a series of steps (Fig. 4). The length of the hypotenuse, 
represented in this way, is 2. As the number of steps N is increased, this length remains 2, even in the limit where 
N — > oo. The fundamental problem is that the hypotenuse is not the N — > oo limit of the zigzag curve. Similarly, it is 
impossible to represent accurately a smooth surface or curve on a cubic lattice. While smoothing increases the length 
scale of curvature for the smooth surface (or, as in Fig. 4, decreases the effective step size), the lattice representation 
of the surface never reaches the continuum limit. 



G 



FIG. 4. An illustration of a classic geometric fallacy: as the number of steps goes to infinity, the length of the hypotenuse 
remains at 2, rather than approaching the correct value of 2 1 ^ 2 . 

Armed with this information, we can now calculate the expected ratio of the continuum value for (A/V) given in 
equation (|l^) and the lattice value in equation ( |i"7j ) . Suppose that we represent a locally flat surface with area A and 
normal vector n on a cubic lattice. The area A then forms a triangular face of a pyramid with sides of length x, y, 
and z. The volume of the pyramid is V = \An = ^xyz, so that A — ^xyz/n. The lattice value for this surface area 
is just Aiat = h(xy + xz + yz). Then (A/V)i a t/ (A/V) is just the mean value of the ratio of Ai at to A, which is 



(A/V)i a t I A 



Mat 



(A/V) \ A 

\{xy + yz + xz) 



\xyz/n 

= (cos a + cos f3 + cos 7) , (18) 

where a, (3, and 7 are the angles between n and the x, y, and z axes. Since n is distributed isotropically, we have 
(cosa) = (cos/3) = (COS7) = 1/2, and (A/V)i a t/(A/V) — 3/2, which is exactly the value obtained by comparing 



equation (fL3j) to equation (17) 



In fact, the solution to this problem for the case of domain walls was noted by Press et al. |lfj and has been used in 
all domain wall simulations JlCj ] - p2 |. In the Press, et al., procedure, the area of a domain wall is simply multiplied 
by the factor l/(\ cosa| + | cos/3| +Tcos7|), resulting in the correct continuum value for the domain wall area. 

A similar weighting procedure cannot be applied in a straightforward way to simulations of cosmic strings and 
monopoles. For strings, a similar calculation yields (L/V)i at /(L/V) — ((x + y + z)/y/x 2 + y 2 + z 2 ) = |, while the 
actual ratio from equations ( |l3| ) and ( pj| ) is (L/V)i at /(L/V) = 6/tt = 1.91, which is slightly larger. For monopoles, 
our geometrical argument fails completely. There is no obvious reason that (N/V) should be any different in a 
continuum or lattice calculation, since the monopoles are geometrical points. Nonetheless, equations ( [l3"|) and jl^ ) 
yield (N/V)i at /(N/V) = 6/n, as in the case of strings. 

This secondary lattice dependence is a result of the way in which we have simulated the string and monopole 
formation. Both the strings and monopoles are treated as arising from intersections of the surfaces defined by the 
zeros of 2 or 3 independent Gaussian fields. Even when smoothed, these surfaces retain a residual jaggedness, as 
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we have already emphasized. The intersection of two such jagged surfaces produces spurious additional intersections 
which would not be present in the continuum limit. As an example, consider the hypotenuse in Fig. 4. When 
represented as a series of steps, two such lines can have a large number of intersections, even when the two lines are 
parallel, while in the continuum limit, there would be at most one intersection, and none in the case where the lines 
were parallel. This secondary effect does not enter into the domain wall calculation, since the domain walls arise from 
the zeros of a single field. 

This problem can be resolved by linearly interpolating the values of each field on the lattice, effectively representing 
the zeros of the field as polyhedral, rather than cubic, surfaces. We have tested such a procedure for the case 
of monopoles. We derive an equation for the value of <j> for each of the three fields in each cubic cell by linearly 
interpolating on the values of each field at the four points (0,0,0), (1,0,0), (0, 1,0), and (0,0, 1). We then place a 
monopole in the cell if the <j> = planes of the three fields intersect inside the cell. The results are shown in Fig. 5. 
They agree extremely well with the continuum predictions, despite the fact that our linear interpolation is still an 
approximation to the true continuum field. 




FIG. 5. As in Fig. 3, but with linear interpolation on the three fields used to determine the monopole positions. 

This procedure can also be used for cosmic strings. In this case, we take the values of our two fields to lie on the 
vertices (rather than inside the cells) of the cubic lattice. For each face of the cubic lattice, and for each of the two 
fields, we use bilinear interpolation on the four values of the field bounding that face to calculate the location of 
the zeros of the field on that face. This gives the intersection (if any) of the two surfaces (f>i = and <f>2 = with 
the given lattice face. If both surfaces do intersect a given lattice face, then we obtain two curves (one for each field) 
which give the intersection of the 4>i = and 02 = surfaces of that field with the lattice face. The intersection of 
these two curves gives the point at which a string passes through that face of the cubic lattice. We then connect these 
intersection points with straight segments (now no longer parallel to the lattice edges) to obtain the string network. 
(The general expression for the location of a single = curve on a lattice face, obtained from bilinear interpolation, 
is axy + bx + cy + d = 0, where a — d are constants and the face is taken to lie in the x — y plane. Hence, this procedure 
can produce two intersections of the <f>i = and 4>2 — curves on a single lattice face, corresponding to two strings 
passing through a single face of the lattice. In this case we treat the string as a single small loop with length equal to 
twice the distance between the intersection points. Such tiny loops become relatively less important as the smoothing 
length increases). The resulting values for L/V as a function of smoothing length are displayed in Fig. 6. Again, 
agreement with our continuum predictions is excellent. 
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FIG. 6. As in Fig. 2, but with linear interpolation on the two fields used to determine the cosmic string positions. 

We can use this linearly- interpolated string network to calculate a "lattice- free" estimate of /oo, the fraction of 
string length in the form of infinite strings, which we define, as in Ref. to be strings which cross the entire 
simulation volume in either the x, y, or z direction. We have calculated foe as a function of smoothing length for our 
linearly-interpolated string simulations. The results are given in Fig. 7. 
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FIG. 7. The fraction of total string length in the form of infinite strings, /oo, as a function of smoothing length, ro, for 
Gaussian smoothing, where linear interpolation on the two fields is used to determine the cosmic string positions. 

To obtain an estimate of the lattice- free value of /oo, we want to choose values of ro which are large enough 
compared to the lattice spacing to eliminate lattice effects, but small enough compared to the total box size to avoid 
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large fluctuations from one run to the next. The results shown in Fig. 6 suggest that the range 3 < tq < 7.5 should 
be suitable. When we perform a weighted average of the values of /oo in this range, we obtain: /oo = 0.716 ± 0.015. 
This result represents a true lattice-free estimate of , but we do expect it to depend on the actual window function 
used for smoothing. Our lattice-free results are strikingly similar to earlier results obtained with Gaussian smoothing 
in a simulation where strings were constrained to lie on a cubic lattice: /oo = 0.71 ± 0.01 ||. It would appear that 
/oo, unlike (L/V), is not much affected by the use of a cubic lattice. 

For the case of monopoles, one interpretation of our results is that our simple-minded lattice simulation leads to 
the formation of spurious closely-spaced monopole-antimonopole pairs. To test this hypothesis, we must distinguish 
between monopoles and antimonopoles in the simulation. Consider first the case where the <pi = surface lies in the 
y — z plane, the </> 2 = surface lies in the x — z plane, and the (pa = surface lies in the x — y plane, so that the 
only non-zero derivatives of the field are 4>i tX1 4>2.y an d 4>3. z - in this case, we define zeros of the three fields to be 
monopoles if 4>i,x4>2,y4>'3,z > 0, and antimonopoles if this product is negative. For other field configurations, the sign 
of the product of the derivatives which corresponds to a monopole changes each time we alter the direction of one 
pair of zeros. For instance, if the cf>x = surface lies in the x — z plane and the 4>i = surface likes in the y — z plane, 
but the 4>s = surface remains in the x — z plane, then 4>i, y 4>2,x4>3-z is now negative for monopoles and positive for 
antimonopoles. Similarly, if we now permute the directions of the zeros of the <p2 and <p3 fields, then 4>i l p<p2,z4'3,x is 
positive for monopoles and negative for antimonopoles. 

Using this definition of monopoles and antimonopoles, we have checked the consequences of allowing closely- 
separated monopole-antimonopole pairs in our simulation to annihilate. We use the same monopole simulation 
as in Fig. 3, but now for each monopole in the simulation, we check for the existence of any antimonopoles 
within an "annihilation distance" ta- If an antimonopole is found within this distance, both the monopole and 
antimonopole are removed from the simulation. Thus, in our final configuration, all remaining monopoles are sepa- 
rated from antimonopoles by a minimum distance of ta ■ We have fixed the smoothing length tq and calculated the 
monopole+antimonopole number density (N/V) as a function of the annihilation distance ta- The results are shown 
in Figs. 8 and 9 for smoothing lengths of 6 and 8 lattice spacings, respectively. The monopole density (N/V)rQ 
decreases as a function of ta, and it must obviously go to zero as ta reaches the size of the simulation volume. 
However, it appears that as ta approaches ro, a plateau value for (N/V)rQ is reached. Moreover, the plateau value of 
{N/V)tq lies at the continuum prediction of equation (|l3|). The statistics are rather poor, but our results do suggest 
that one interpetation of the spurious monopoles in the lattice simulations is that such monopoles are produced as 
closely-spaced monopole-antimonopole pairs. 
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FIG. 8. The total number of monopoles+antimonopoles per unit volume, (N/V), multiplied by the cube of the smoothing 
length ro, as a function of the "annihilation distance" ta, where monopole-antimonopole pairs separated by a distance less than 
ta have been removed from the simulation, and ro is fixed at 6 lattice spacings. The solid line gives the analytic prediction for 
(N/V)rQ in the continuum limit. The dotted line is the lattice-based analytic prediction. 
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FIG. 9. As in Fig. 8, with ro = 8 lattice spacings. 



IV. DISCUSSION 



Our results indicate that lattice simulations for the formation of topological defects must be used with some caution. 
Even when the Gaussian fields used in these simulations are smoothed to reduce lattice effects on small scales, a residual 
lattice dependence remains. Another way to express this result is to note that smoothing eliminates the small-scale 
inhomogeneity introduced by using a lattice, but it cannot eliminate the anisotropy. Strings, for example, simulated 
on a cubic lattice are required to lie only in the x, y, or z directions, while a genuine continuum simulation would 
allow an arbitrary direction for the string segments. 

For the case of domain walls, a simple multiplicative factor which depends on the direction of the domain wall 
can be used to eliminate this problem; all domain wall simulations have used such a factor |l(| - |l2| . For strings, 
the correct multiplicative factor does not amount to simply averaging over the directions for a straight string: in the 
lattice simulations there is residual small-scale structure along the strings and probably an accompanying distribution 
of tiny loops (of order the lattice spacing). These effects change the multiplicative factor from 1.5 (its value when 
only anisotropy is taken into account) to 1.91, but they should not affect the length distribution of loops (except for 
the smallest loops). In practice, nearly all previous studies of cosmic strings have used a simple lattice simulation of 
the kind described here Q - Q , as have some monopole simulations |l3| (although see Ref . [[TiJ for a continuum field 
treatment of monopoles) . 

For cosmic strings and monopoles, the problem can be solved by linearly interpolating the values of the field on 
the lattice. This procedure yields the correct continuum value for the monopole and string densities; for the case 
of monopoles, a similar result can be achieved in a simple-minded lattice simulation by eliminating close pairs of 
monopoles. For the case of cosmic strings, this interpolation scheme produces a truly lattice-free estimate of the total 
string length in the form of infinite strings: = 0.716 ± 0.015. This is nearly identical to the value obtained from 
a simple lattice simulation (with Gaussian smoothing) where the strings were constrained to lie on the edges of the 
lattice [^| . Apparently foo is much less sensitive to lattice effects than is the total string density. This is not surprising, 
since the effect of the cubic lattice is to multiply the string density by a geometric factor; if this factor is the same 



11 



for both closed loops and infinite strings, then /oo is unaffected. However, even in these lattice-free simulations, we 
expect the value of /oo to depend on the window function used for smoothing (as in reference || ) . 

This linear interpolation procedure can also be extended to domain walls, providing an alternative to the weighting 
factor used in references Jl(J - jl2| . In practice, these two methods should give nearly identical results for any domain 
wall quantities of interest. 
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